NTUA-92/00 



QED in a Strong External Magnetic Field: 
Beyond the Constant Mass Approximation 

J.Alexandre, K. Farakos, G. Koutsoumbas 
Department of Physics, National Technical University of Athens, 
Zografou Campus, 157 80 Athens, GREECE 



Abstract 

We solve the Schwinger-Dyson equations for QED in 2+1 or 3+1 dimensions in the 
presence of a strong homogeneous external magnetic field. The magnetic field is assumed 
strong enough, so that the lowest Landau level approximation holds, but the usual assumption 
of a momentum-independent self-energy is not made. In 2+1 dimensions, the scaling with 
logarithm changes to a square root dependence on the magnetic field, but the most spectacular 
result takes place in 3+1 dimensions, where the constant mass approximation turns out to 
be unreliable and the (momentum-dependent) dynamical mass is larger by several orders of 
magnitude compared to what has been found till now using the constant mass approximation. 
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1 Introduction 



Magnetic catalysis of chiral symmetry breaking is well known by now as a concrete example of 
dynamical symmetry breaking, both in 2+1 and 3+1 dimensions. It is, of course, interesting 
from the field-theoretical point of view to have an example of dynamical symmetry breaking 
where explicit calculations may be performed. The non-perturbative tool used has been a 
consistent truncation of the SD equations based on the domination of the dynamics by the 
lowest Landau level (LLL) in the regime of strong magnetic fields. The results have shown that 
the presence of a homogeneous external magnetic field is enough to trigger dynamical mass 
generation even at very weak attractive coupling between the fermions. 

The long history of the subject involves among others the very interesting issue of its 
role in models of the QCD vacuum |jl| and a possible involvement in a new phase of QED with 
spontaneous chiral symmetry breaking, which has been claimed to offer an explanation of the 
e'^e~ pairs observed in the heavy ion experiments in GSI 0. 

The Standard Electroweak Model in the presence of a constant magnetic field has been 
examined by Ambjorn and Olesen 0. Their result is very briefly that, for strong enough 
magnetic field strength, they obtained a W and Z condensate solution and a lattice of vortex 
lines is formed, while at even stronger magnetic fields a phase transition to the symmetric 
phase takes place; the critical temperature needed for the symmetry restoration is substantially 
smaller than the one without magnetic field. Thus, it is seen that the magnetic field may induce 
a transition from the broken to the symmetric phase of the Electroweak Standard Model. These 
ideas of symmetry restoration in the presence of a magnetic field are discussed also earlier (see 

!)• 

The ground state of the electroweak theory has also been studied in the presence of a 
hypermagnetic field to the end of stregthening the electroweak phase transition so that it might 
give rise to baryogenesis within the Standard Model |Q. Moreover, the numerical work done 
and the contribution of ring diagrams in the effective potential show that the hypermagnetic 
field strengthening of the phase transition is not enough to satisfy all the standard model 
baryogenesis conditions. 

Models for the formation of the strong magnetic fields that have been observed on the 
scale of galaxies are strongly based on the existence of extremely intense primordial magnetic 
fields [^]. These astrophysical and cosmological considerations are another field where magnetic 
field induced chiral symmetry breaking is relevant. 

The model of a QED field theory with fermions has been suggested to describe the long- 
wavelength limit of the high-T^ superconducting cuprates [0]. Aspects of the behaviour of these 
materials, such as the scaling of their thermal conductivity, may find their explanation in the 
framework of the model under the influence of an external magnetic field 

The theory of dimensional reduction and magnetic catalysis has been developed by 



several groups |T0| [IT], |12[. They have treated the problem of a constant external magnetic 
field acting on either Nambu Jona Lasinio model or QED in 2+1 or 3+1 dimensions. It 
was found that the magnetic field is actually a strong catalyst and even the weakest possible 
attraction between the fermions is enough for dynamical mass generation. In addition, in 2+1 
dimensions, mass is generated dynamically for any number of flavours and not only for < 4, 



as is the case without magnetic field and the critical temperature for symmetry restoration 
is much smaller than VeB. 



There have also been found solutions [|1J] of the Schwinger- Dyson equations in the limit 
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of weak magnetic field for (2+l)-dimensional QED with both quenched and dynamical fermions 



14| . Also, if there is a coupling of the fermions to scalar fields in the presence of a magnetic 
field, a strong fermion condensate appears analogue to a/ eB ||15|| . 

All of the above approaches have used the constant mass approximation, that is the self 
energy was supposed to be independent from the momentum. In this paper we try to make a 
step towards discarding this rather crude assumption and see the consequences in both the 2+1 
and the 3+1 dimensional theories. To be able to proceed, we restrict ourselves to the strong 
field regime. 



2 Fermions in a constant magnetic field 

To fix our notations we shortly review here the characteristics of fermions in a constant external 
magnetic field, both in 2+1 and 3+1 dimensions. 

The model we are going to consider is described by the Lagrangian density: 

jO- = —U.F" + t^D^Y"^ - m^^, (1) 

where = + iga^ + ieA^^^, is an abehan quantum gauge field, f^i, is the corresponding 
field strength, and A^^^ describes an external electromagnetic field; in this work it will be a 
homogeneous magnetic field, constant in time. Notice that the fermions feel both the quantum 
and the external gauge fields, however we have allowed for different coupling constants, g and 
e. In three dimensions it makes sense to keep different couplings 0, while in 3+1 dimensions 
g is set equal to e. We recall the usual definition = Ana. 
We will choose the "symmetric" gauge 

Ar\x) = -^X2, = +|xi, (2) 

with the remaining components of the potential vanishing both in 2+1 and 3+1 dimensions. 
We notice for further reference that this potential has the property: 

xf'Al'^'iy) + y^Al'^'ix) = (3) 
We know from the work of Schwinger [|T6[ that the fermion propagator is given by: 

S{x,y) =e''''"^^-^'^y^S{x-y), (4) 
where the translational invariant propagator S has the following Fourier transform: 



S{p) 



dsexpis — Pj_ 



tan( ei? 






eB 


s 



m 



(p"7" + m){l + Yt ts.n{\eB\s)) - p_l7_l(1 + tan^(|e5|s)) 



(5) 



where p-^ = p" = p° in 2+1 dimensions and = {p'^,p^) in 3+1 dimensional space 

time. Similar notations hold for the 7-matrices. 

There exists an additional representation for the fermion propagator, namely the Landau level 
expansion |jl^. This expansion is particularly useful in the strong field regime, \eB\ » rn^, 
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since one may keep only the first term and get the so-called lowest Landau level approximation 
(LLL) 



SLLLip) = .e-^i/l^^iq^ (l - t^'^hgieB)) (6) 
It is this strong field regime that we will consider in this paper. 

3 Schwinger-Dyson equation with an external magnetic 
field 

The Schwinger-Dyson equation for the fermion propagator in d dimensions is given by the 
expression: 

G{x,y) = S{x,y) - Aira J d'^zid'^Z2d'^Z3d'^ZiS{x, Zi)-f''G{zi, Z2)T''{z2, Z3, Z4)G{zi,y)A^^{zs, Zi), 

(7) 

where G and A^j, are the full fermion and photon propagators and the full vertex. We show 
now perturbatively that G factors in the same way as the bare propagator S in equation (Q), 
in other words it can be written in the form: 

G{x, y) = e'^"'^M^*(^)G'(x - y). (8) 

Making a loop expansion in (|^ and keeping only the tree level and the one loop terms we can 
write: 

G{x, y) = S{x, y) - 47r« j d^z^d^z^Six, zi)-i^S{zi, 22)7'^(^2, y)D^,{z2, z^) + ... (9) 

Note that we have used the bare vertex T^{z2, 2:3, Z4) = 7^5(^2 — Z3)6{z3 — Z4), and that to this 
order in the loop expansion only the free photon D^j, and fermion S propagators appear. 
Let us now isolate the one- loop term in the expansion (Qj. This term reads: 

G^'\x,y) = -Alia j d''zid''z2S{x,zi)^>'S{z^,Z2)YS{z2.y)D,,{z2,z^) 

= -Ana I d''zid''z2~S{x - ;2i)7^5(2i - ;22)7"^(^2 - y)D^,{z2 - z{) 

X expie + A;^^*(22) + zi^A^^\y)) (10) 

Using the change of variables Zi = y — Z2 = x + v and writing r = x — y, we obtain: 



G(i)(a;, y) = -Ana J d'^ud'^vS{r + u)^''S{-r -u- v)YSir + v)D^^{r + u + v) 

X exp ze - u) + (y - uYAfix + + (x + vfAfiy)) (11) 

Taking advantage of the property (^) of the potential A^^*, we finally obtain: 
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G^^\x,y) = e^^"''^r(2/)G(i)(r), (12) 

where 

G(i)(r) = -Ana J dVve'^"'^'^^'^''^ S{r + u)YS{-r -u- v)YS{r + v)D^^{r + u + v) (13) 

which is exactly equation (||). To be precise, the full propagator equals a phase times a trans- 
lational invariant quantity. Proceeding the same way we can see that all the higher loop 
corrections behave as the one loop, showing perturbatively that G factorizes in the same way 
as S. 

We can also see that the fermion self energy S = G^^ — S^^ has also the same structure 
as S", since the Schwinger-Dyson equations can be written as: 

E(x, y) = 4^7^^ j d^Zid^Z2G{x, ^1)^(2/, Zi, ^2) A^,(z2, x), (14) 
such that the one-loop self energy is: 

EW(x,y) = ATiaYG{x,y)-i''D^,{y,x) 

= e''=""^^"(^)E(a;-2/) (15) 

where S(r) = ATTa'y^G{r)'y'^D^iy{r). The higher order corrections to S will have the same phase. 

Since we will deal with the LLL approximation only, we cannot define the inverse of the 
fermion propagator (|^) because of the projection operator (1 — 27^7^)/2; thus we cannot define 
the LLL approximation for S(x,?/). For this reason we will prefer to use the Schwinger-Dyson 
equation (|^) and not (p!^. We will also take the bare vertex, so that our starting equation is: 

G{x, y) = S{x, y) - A-aa f d'^zid'^Z2S{x, zi)j''G{zi, Z2)YG{z2, y)A^,{z2, zi) (16) 



We rewrite equation (|T6D in terms of the translational invariant fermion propagators S{r) and 
G{r) 

G{r) = S{r) - 47ra j d''zid''z2S{r - zi)YG{zi - Z2)YG{z2)A^,{z2 - ^i)e^-('-^2)"^r(-i) (17) 
where r = x — y. Making a Fourier transform and performing the integration over Zi we obtain: 

G{k) = S{k) -Anai{k), (18) 

where 

m = jd'zj ^^-^^~S{p-q + eA^^\z))^^G{p)^^G{k-eA^^^ (19) 

We will look for a solution of ( [T^ ) where the translational part of the full propagator reads in 
the spirit of the LLL approximation (equation (^): 
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G,Up) = .e-^i/l-^l ^g^^_+^^ (l - ^JVs9ieB)) , (20) 
where Zp is the wave function renormahzation and Mp the dynamically generated mass. 

4 Magnetic catalysis in 2+1 dimensions 
4.1 Integral equations 



We will look for a solution of (p!8| ) with the expressions (g) with m = and ( pO| ) substituted for 
the free and the full fermion propagator respectively. We are using four-component representa- 
tion for the fermions and the 7 matrices. To obtain the equations that Z^ and must satisfy, 
we multiply I{k) respectively by {Zkk'^'y^ — Mk) and 7°(Zfc/c''7° — Mk). We do a Wick rotation 
and then perform the integration over zq = iz^, which leads to a factor 27i6{q3 — ps + k^). We 
obtain then the equations: 

where = iz±{q±-p±+k_i_)-j^^ [{p± - q± + eA''''\z)Y +p\ + {k± - eA^^^z))"^] and T>{qs, q±) = 
—iDoQ^qE). We note that only the component Dqq of the photon propagator plays a role in the 
Schwinger-Dyson equations, due to the projection operator (1 — z7^7^)/2 which appears in the 
fermion propagators when the LLL approximation is used. 

The integrations over pj_ and z± are straightforward. We also have the auxiliary relations: 



tr 



{{Z.zksi' - M,) (G{k) - S{k))} = 4ze-'=i/l^^l(l- 
tr {J'{z,^ksJ'-M,){G{k)-S{k))} = 4e-^i/l^^l^ 

such that (|T^) yields the final integral equations: 



.{l-Z^) = - dv due~-''-^^V{v-.,u) (21) 

TT J -00 JO ' H'v 



We have introduced the dimensionless variables a = a/ J\eB\, u = Jq\/\eB\, v = p^/ J\eB\, 



K = fcs/y |e-B| and /i^ = M^/ y\eB\. The reader can check that the integral equations ( pT| p^ ) 
are consistent with the ones in [|12| . We note that /i^ and Z^, depend only on the dimensionless 
ratio a. 
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The free photon propagator A^,^ will be used, since we show in the appendix that it 
receives no corrections in this case. In particular, it is shown that the one-loop polarization 
tensor is suppressed by a\eB\~^^'^, such that in the LLL approximation we will use in (^, ^) 
the Euclidean photon propagator 

V(V -K,U) = — TTT 1 - ^r^— — '—^ + A ^ ^ 23 

where A is the gauge fixing parameter. 



4.2 Constant mass approximation 

Before proceeding with a detailed analysis of ( PT| , we start with a first approximation 
which consists in setting = 1 and a constant dynamical mass /Xc in the integral equations 
so that in the Feynman gauge A = 1 equation (22) reads: 



const ^'^ r ^ f°° ^^e""'/^ i f24) 

The integration over v is done by the residue theorem and leads to: 

const ^~ r ^^^-uV2 ^ + f^c ^25) 

Thus the constant dynamical mass is obtained by solving the equation = which reads: 



OO " - 1 



rfMe-"/2_^ (26) 

U + He 



The study of ( PB| ) has also been done in 

Putting = 1 as a first approximation and taking a constant mass in equation (^) leads to 
the following second approximation to : 



const 



+ /U 

such that its value for zero momentum is given by: 



We plot the dimensionless dynamical mass gap a ^[i^Zc versus ^\eB\la in figure ^ and com- 
pare with the mass gap obtained with /Uk=o and Z^=<^ solutions of (0,|^) that we find in the 
next section. We can see that the constant mass approximation leads to an over-evaluation of 
the dynamical mass gap and that it gets worse as the magnetic field increases. 
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4.3 Momentum dependent solutions 

We now concentrate on the momentum dependence of the dynamical mass and the wave func- 
tion renormalization to obtain a more precise analysis of the solution of ( PT] , P^ . 
Let us first study the wave function renormalization Z^. We will actually consider the one- loop 
approximation, consistently with our treatment of the Schwinger- Dyson equation (|^), where 
we considered the bare vertex. Since there is no pertubative solution of the integral equations 
for the dynamical mass /x^, our numerical solution of the system of integral equations will cor- 
respond to some resummation of diagrams. The solution will be found through an iterative 



procedure: a trial function will be put in equations (|21| , P2D to yield a first estimate of and 
/ifj- We will then take this output and substitute it back in the right hand sides of the integral 
equations to obtain better approximations. As we will see, this iterative procedure does finally 
converge to a solution. 

The trial function will be given by the solution of a differential equation verified by [i^, that 
will be derived now. This differential equation is derived also in for the (3+l)-dimensional 
model in what they called linearized approximation. 



Let us now proceed with equation (^2j). We start by splitting the integration over v in 
(^) in two parts: 

TT J\v\<K Jo Z^V^ + fif, 

Notice the approximations used: V{v — k,u) ~ V{k,u) for |t>| < k and V{v — k,u) ~ V{v,u) 
for |f I > K. Then the derivative with respect to k, leads to the following equality: 

where we defined: 

^ A ||^°°rf^e-"'/2MP(/€,n)} (31) 
A new derivative with respect to k finally gives: 

/X - />L - -f. , = (32) 



From now on we take Z,^ = 1 in (p^). We should stress that this differential equation is only 
approximate, and the approximation (^) is not controled by any small parameter. Since we 
have no full control on the assumptions made, we consider the solution of (|3^ ) only as a trial 
function which could help us with the solution of the integral equations. It turns out that 
the iterative procedure described above converges quickly yielding a solution for Z^ and /i^ 
reproducing itself after a small number of steps. This may be seen in figure ^ for /i^. We find 
that the convergence of is even quicker. It turns out that the solution found as the limit of 
the iterative procedure does not depend on the trial function that one starts with; however, 
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Figure 1: versus k for a = .03 (d=2+l) 



taking the solution of the differential equation gives a much better convergence, that is the 
speed of convergence is increased. As can be seen on figure |I|, is almost flat for k < /iK=o, as 
was also emphasized in |^ for 3+1 dimensions. 

To solve (|3^), we needed the initial conditions that the function /x^j has to satisfy. Sup- 
pose that /iK=o = /io- In the Feynman gauge, the function has the following behaviour around 

K = 

1 

~ as K (33) 



K 



Thus if we wish /i" to be finite in k = 0, 



implies that: 
2a 



TT/io 



(34) 



Next we have to determine /io- We invoke the fact that /i^ ^ as k — oo and thus choose /iq 
such that limK^oo = 0. As a side remark we note here that the trial function we obtain is not 
an even function of k; this is a consequence of the approximations made to get the differential 
equation (p2|). The even function fi^ is found after plugging the trial function into the integral 



equations (|2T 



I). Equation (|3^) has been solved numerically through the use of a fourth order 
Runge-Kutta method. The estimate of the function as well as its derivative has been done 
using a Gauss-Hermite quadrature of order 40. 

In figure ^we plot the dimensionless mass gap /^k/^k as a function of the dimensionless 
momentum k. One of the curves comes from the solution of the integral equations (using the 
iterative procedure), while for the sake of comparison we also plot the solution of the differential 
equation (0). We note that in the differential equation we have set = 1; on the contrary, 
changes during the iterative procedure for the solution of the integral equations. We observe 
that the quantity fin/Z^ derived from the differential equation is very close to its final value that 
is found after substituting it the integral equations. This explains why the iterative procedure 
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converges fast if we use /i^ as input. It is obvious from the figure that the mass gap decreases 
rather rapidly with increasing momentum and takes its maximum value for n = 0. Thus it is 
to be expected that we will have big differences from the constant mass approximation, which 
assumes the same value for < k < oo. The small momenta will be given the major role, in 
accordance with the fact that mass generation is an infrared phenomenon. 

In figure ^ we try to get contact with |^ and plot the dimensionless dynamical mass 



ji/a = M/a versus \J\eB\/a for two cases: the constant mass approximation obtained from 
(p6|) and the momentum dependent version, coming from the solutions of the integral equations 



OH)- 

In the constant mass case, we find a logarithmic fit, which reads: 



fit V 

a ficonst = n^o 111 with uiq ~ .765 ± .001 (35) 



a 

or 



donst = mo ^ — In ^ '- (36) 

\eB\ " 



a 



eB\ 



in agreement with the result of [§. 

For the momentum dependent solution, we choose to plot the value of this dynamical 
mass at K = 0. versus \eB\. The fit now is dramatically modified and reads: 



1 ft A/ 

a-Vsin = mi+ m2- with mi = 2.65 ± .01 and ma = (2.1 ± 0.1) x 10"^ (37) 

a 

Of course, changing a logarithm to a square root is a fundamental difference from the previous 
case. We also observe the quantitative change in the value of the mass in comparison to the 
constant mass approximation: for the magnetic fields which are plotted in figure ^ the mass 
is smaller by a factor of two to three. We should stress that the above results have been 
produced using the lowest Landau level approximation, which holds provided the magnetic 
field is strong. This means that our results, shown in figure |^, should not be trusted in the 
region near the origin of the horizontal axis. In this respect we observe in the figure that the 
data for the momentum dependent mass are well represented by the above fit in the region 



of large magnetic fields down to the value y\eB\/a ~ 100. We may interprete this change of 
behaviour of the data for small magnetic fields as been due to the LLL approximation, which 
is poor in this region; this may yield an estimate of the limit of its validity. We note here that 
with zero magnetic field, the dynamical mass obtained by the Schwinger-Dyson equation in the 
constant mass approximation is bigger than the mass obtained solving the equations using a 



momentum dependent self energy for the fermion [13 



4.4 Fermion condensate 

The fermion condensate is given by 



0|^^|0) = limtrG(x,y) = ti J ^^G{p), (3 



(2vr) 
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Figure 2: fu,^/Z^ versus k, ioi a = .001 (d=2+l) 
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;ure 3: Dimensionless mass gap fj,/a versus J\eB\/a {d— 
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Figure 4: Dimensionless fermion condensate < 0|'?/''?/'|0 > /\eB\ versus J\eB\/a (d=2+l) 



which, in the LLL approximation, leads to: 



-Ai 
leB 



(2vr)2 



dv 



-pl/\eB\ 



/it. 



idps 



271 zy, + M^ 



2tt'^ J-oo Z^v"^ + nl 

In the constant mass approximation with = 1 the integral (|39D yields the simple result: 



(39) 



\eB\ 



(40) 



It is worth noting that in this approximation the condensate does not depend on the dynamical 
mass The condensate in the constant mass approximation grows linearly with \eB\ with 
a slope equal to l/(27r) ~ 0.1592. We compared (^D|) with the condensate given by equation 
(^) using the momentum dependent functions /i„ and Z^ found numerically in the previous 
subsection. We show in figure I 



the dimensionless condensate {Q\ipip\Qi /\eB\ versus \J\eB\/a, 
both for the momentum dependent and the constant mass solutions. We see that the condensate 
varies linearly over a wide range of 
vanish in the limit \eB\ 



\eB\/a 
in 2+1 dimensions 



We remark that the fermion condensate does not 
but we cannot explore this limit within 
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the LLL approximation, which is only reliable for strong magnetic fields. 



5 Magnetic catalysis in dimension d=3+l 
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5.1 Integral equations 



The techniques used in the (3+l)-dimensional model are similar to the ones developed in the 
previous section, so we will only give the results. 

Starting again from the Schwinger-Dyson equations we obtain in the same way as before 



= -,jdvj^ due ' ^^^V{v-.M 



/X. = -Jd\r due--'l'-P!^V{v-RM (41) 
TT^ J Jo Zf,v^ + af, 



where now k is the two dimensional vector representing two components of the dimensionless 
Euclidean momentum: 

K = }. — (A;4, ks) = }- — {-iko, ks) (42) 
\eB\ J\eB\ 



In 3+1 dimensions, the corrections to the photon propagator are not suppressed by inverse 
powers of the magnetic field, so we use the expression [|]: 

where n(g5_,Q'|) is the polarization tensor and (^Jl^ and g^^, represent the restrictions of the 
metric tensor to the (3, 4) and (1, 2) directions respectively. Because of the projection operator 
(1 — i'~f^'~f'^sg{eB))/2, appearing in the fermion propagator in the LLL approximation, only the 
part of D^i, containing gfjl^^ will contribute to the Schwinger-Dyson equations. Therefore we will 
only have to consider the expression 

where we will use for the polarization tensor the form |Tl[| : 



gjn(gj,gi)^-^e-i/(^l^^l), (45) 

II M ^ 



a being the QED coupling renormalized at the scale y\eB\. Actually the approximation (|g) 
is valid only for <lf » M"^ where M is the dynamical mass of the fermions, but the latter will 
be very small compared to the typical momenta we will consider, so the approximation is good. 
Therefore the photon propagator in dimension 3+1 will be: 

Di^^iq) = -I _ 2a|g^U-gi/(2|eB|)' (^6) 
H -j^ \ I 

such that the quantity V{v — k, u) in (^Tj) will be replaced by: 

^("^ - ^' = «^ + (.--^)'+^e-«V2 (47) 
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The integral equation for /x^ is consistent with the results of |^ where the authors do not consider 
the contribution of the wave function renormalization. It is known that in the dynamical 
symmetry breaking in QED, without external field, is equal to one, if we employ the Landau 



gauge |]T9|, sect. 8. 4]. However, with an external magnetic field, the photon propagator (|46|) is 



such that we find 7^ 1 as can be seen in figure ^ 



5.2 Constant mass approximation 



We first have a look at the constant mass approximation of (0), to facilitate comparison with 
0; we do not consider in this subsection the corrections to the photon propagator. Denoting 
p = K^, we consider the integral equation for p = and Zy = 1 and set fip = fiy = fic- The 
integration over v then leads to the following gap equation: 



a 



TT JO 



u 



We plot in figure ^ the dynamical mass obtained from 
of [^, which is: 




(4J 



as well as the analytical estimate 



analyt 



C exp 



TT / TT \ 
'2 



(49) 



and choose C = 1. We can see that this estimate agrees very well with our results. But we 
will see in the next paragraph that the momentum dependent dynamical mass gives completely 
different results from the constant mass approximation in 3+1 dimensions, which was not the 
case in 2+1 dimensions. From the numerical point of view, this comes from the fact that in 
3+1 dimensions the dimensionless dynamical mass pc is not present in the left-hand side of 
the gap equation 



contrary to what happened in (|26|) , the corresponding equation in 2+1 
dimensions. The result is that the order of magnitude of the integral should be 1/a in 3+1 
dimensions. This pushes the integral to very big values, that can be obtained by extremely 
small values for p^. On the other hand, such problems are not present for the momentum 
dependent solution; this leads to a serious discrepancy, by several orders of magnitude, between 
the momentum dependent solution of (|41|) and the constant mass approximation which therefore 
is not reliable in 3+1 dimensions. We note that this fact is independent of our taking into 
account Z^, since the latter is very close to 1 and thus cannot give a change of several orders 
of magnitude (see figure We also note that such problems do not arise in 2+1 dimensions, 
since, by equation (p6D , the order of magnitude of the integral is Pc/a in this case, so the 
smallness of a is compensated by a small number, such as pc; thus moving to a momentum 
dependent solution in 2+1 dimensions brings in some quantitative differences but not the huge 
ones characterising the 3+1 dimensions. It seems that the constant mass approximation in 2+1 
dimensions is relatively reliable. 



5.3 Momentum dependent solutions 

We are looking for solutions of (^Ij) which are even functions of the momentum. Thus we 
consider p and Z as functions of the variable p = = k'^/\eB\. The differential equation 
similar to ( p^ ) is then obtained in the same way and reads: 
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Figure 5: Dimensionless dynamical mass fic versus a in the constant mass approximation 
(d=3+l) 



a 



where the function /„ is given by 



d 

dp uo 

Equation ( |5DD is consistent with the one given in 



roo „ 

due-"" /\V{^,t 

and we will set Zn 



(50) 



(51) 



1 to solve it. 



Coming back to the integral equations (^T]), we can perform the w-angular integration if we 
look for functions Pp and Zp depending on p only. The result reads: 



a 



271 p Jo 



a f°° , 
— / dv 

TT Jo 



dv / due-^ ^'^^ 
Jo V + pf, 

upv 



Z^v + pi 



+ P + V + 2a/7rexp(-uV2) 

(m^ + P + V + 2a;/7rexp(— ^2/2))^ — Apv 
1 



, (52) 
[v? + p + V + 2a/TT ex.Y>{—v? /2)) —Apv 

where we have used the one-loop approximation Z^^^ for the wave function renormalization Z^. 
We solve numerically (|S^) in the same way as we solved the equivalent set of equations 



in dimension 2+1. We first solve the differential equation (|50D and use the solution as a trial 
function for the integral equations ([52|) . Then the outcome is used as feedback in the integral 
equations (^ and the resulting iterative procedure is converging, somehow slower as compared 
to 2+1 dimensions. To find the initial conditions of pp satisfying (^0|), we note that: 



1 

Yp 



when p — > 0, 



(53) 
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Figure 6: Wave function renormalization Zp versus p for a = .01 (d=3+l^ 



so that (|50|) implies that: 



Of 



/i 



p=0 



2vr/io 



with 



(54) 



if we demand that fip is twice differentiable as p ^ 0. We then find Hq by imposing the condition 
hmp_,oo/ip = 0, as we did in 2+1 dimensions. 

In figure ^ we depict the evolution of the dimensionless dynamical mass at zero momentum, 
/iK=0) with the coupling a, as well as a fit which is given by: 



..fit 

pp=0 



with /io = .0300 ± 0.0001 (55) 

The dimensionless mass gap {pp/ Zp)f^=Q follows the same fit, since ~ 1. We may compare 
figures | and |^ and see the huge difference (by some orders of magnitude) between the constant 
mass approximation and the momentum dependent solution in 3+1 dimensions. The first 
conclusion is of course the that one cannot really trust the constant mass approximation, as 
already stated above. However, such a change will make itself felt and, if this analysis is 
correct, the relevant chiral symmetry breaking triggered by external magnetic fields should 
have measurable physical consequences. 



5.4 Fermion condensate 

The fermion condensate in the LLL approximation is given by: 

Since we compute the condensate with the momentum dependent mass and wave function 
renormalization, we do not have to compensate any divergence: the condensate is convergent. 
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Figure 7: Dimensionless dynamical mass /ip=o versus a (d=3+l) 



In a perturbative expansion, the condensate has to be renormahzed to any order but here 
we started from the Schwinger- Dyson equation which is non-perturbative and leads us to a 
convergent resummation of graphs for the condensate. We can note that the divergence in the 
computation of the fermion condensate for the constant mass approximation can be removed 
if we substract the condensate for \eB\ = 0. We obtain then a convergent expression which 
vanishes as — 0. To see this, let us go back to the Schwinger representation of the 
fermion propagator and take away the gauge dynamics, putting a 0. The computation of 
the condensate is straightforward and leads to: 



OlibihlO 



OI^^IO 



m\eB\ r°° ds 



|e_B|=0 



47r2 



coth s 



\eB\ 



m 



s\eB\ 



(57) 



In our numerical study with zero bare mass, we also find a vanishing condensate for a 0, 
at least for strong magnetic fields since we only considered the LLL approximation. Thus in 
this approximation we do not find any critical coupling for the condensate to be generated. We 
plot in figure ^ the dimensionless fermion condensate Cf = _B~^/^ (olipiplOj versus the coupling 
a, taking into account the relation = Awa. 

We note that the analytical approximation of the fermion condensate given in [^] does 
not fit to our numerical results. The difference lies in their taking the constant mass approxi- 
mation and a cut-off \eB\ in the momentum space. Instead we found the analytical fit: 



Jit 



Co a 



1/4 



with Co = .00590 ± 0.00001 



(5^ 



for the dimensionless fermion condensate. 
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Figure 8: Dimensionless fermion condensate B < OlV^-i/^IO > versus a (d=3+l) 



6 Conclusions 

Our results for the (2+l)-dimensional model are relevant for condensed matter Physics, in 
particular models for high-Tc superconductivity, while the (3+l)-dimensional case has to do 
with the Electroweak Phase Transition in the Early Universe. 

• As we already stated in the beginning, there has been a suggestion |^ that the high-T^ 
superconducting materials can be described by an effective theory, namely a relativistic 
gauge field theory in 2+1 dimensions. The complete model was SU{2) x f/(l) symmetric 
with a doublet of massless fermions p2|. However, simpler versions, such as QED in 
2+1 dimensions have been studied in the same context. The theory describes the quasi- 
particle excitations about the nodes of a d-wave superconducting gap. The Fermi velocity 



coincides with the effective velocity of light. In a recent experiment p3| it has been found 
that for strong magnetic fields (about 0(1) to 0(10) Tesla) there appear plateaus in the 
thermal conductivity of the material, signalling the opening of a new gap at the nodes, 
induced by the magnetic field. The critical temperature Tcru for the appearance of this 



superconducting gap scales with the magnetic field as Tcru ~ y Assuming that the 
critical temperature (where the dynamical mass vanishes) is roughly proportional to the 
dynamical mass at zero temperature, we see that our equation (|37|), which says that the 



dynamical mass scales as y\eB\, agrees with this experimental result. Of course, a true 
finite temperature calculation must be done if a better founded assertion can be made. 

Concerning the dynamical mass generation and the fermionic condensate in 3+1 di- 
mensions, we mention that around the Electroweak era magnetic fields of the order of 
10^^ — 10^^ are naturally obtained These correspond to a dynamical mass of the order 
of IGeV, if we set 1^ ~ independently from the vacuum expectation value of the 
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scalar field. This is already a big value for the mass with eventually observable contri- 
bution in the effective potential at high temperature. However, if we want to give an 
explanation for the galactic magnetic fields (of the order of 10~^ Gauss), the estimate for 
the magnetic fileds during the Electroweak Phase Transition raises to 10^^ — 10^^ Gauss, 
according to 0. If this is the case, the dynamical mass is of the order of 100 GeV and 
we believe that this possibility deserves further study in connection with Cosmology and 
the Electroweak Phase Transition. 
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A Polarization tensor in magnetic field (d=2+l) 

The free photon propagator in the presence of a magnetic field is ||12| 



DM - [9,. - ^ j + ? [^^■' ^^)^ 1^ ^ ^ 

where A is the gauge fixing parameter. Because of the projection operator (1 — i'y^'y^sg{eB))/2 
that we have in the fermion propagators (§) and (^), the components of D^^, which include 
only the transverse coordinates (1, 2) will not play any role in the Schwinger-Dyson equation. 
Therefore in the LLL approximation the free photon propagator will be 

fli."fe) = -^(i-^)-"^ (60) 

We show now that the corrections to the photon propagator in dimension 2+1 need not be taken 
into account in the LLL approximation. In Euclidean space-time the one-loop level polarization 
tensor is given by p4| , pT 



n^.b) = {p\u - P,Pu)No{p) + (pi^i - pipi)N,{p) (61) 

where 



Noip) = 


a f°° ds 


Jo s 7-1 


Ni{p) = 


a f°° ds 


2y/TT Jo S J -I 


eB\s and 





dve~''^'^° — — — [cosh(2;t>) — t>coth(z) sinh(zi;)] 
sinh(z) 

2z 

dve-"'"' 3— [cosh(2) - cosh{zv)] - Nq {p) (62) 

sinh (z) 



9 1 — v'^ r, coshf^;) — cosh(2;f ) ^ , , 

bo = m' + ——pI + \ ■ . pI 63 

4 2zsmh[z) 
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We will need only the component noo(p) in the LLL approximation, thus we only need to look 
at Nq{p) in the strong field limit \eB\ » s^^. Making an expansion for z >> 1 we obtain 



after the interation over v : 



Jo s ^ 



which finally leads to: 

Noip) 



2a 



irleBl 



Jo 



l + O 



\eB\ 



(64) 



(65) 



We see then that the one- loop polarization tensor vanishes as a|e-B| 
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